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ABSTRACT: 

Based on Refs 1 3,4 1, the unified variable-domain variational theory of hybrid 
problems for rotor-flow 11,2] is extended to fully 3-D transonic rotor-flow with 
shocks, unifying and generalizing the direct and inverse problems. Three 
variational principle (VP) families have been established. All unknown boundaries 
and flow discontinuities (such as shocks, free trailing vortex sheets) are 
successfully handled via functional variations with variable domain, converting 
almost all boundary and interface conditions, including the Rankine-Hugomot shock 
relations, into natural ones. Tins theory provides a scries of novel ways for blade 
design or modif ication and a rigorous theoretical basis for finite clement 
applications and also constitutes an important part of the optimal design theory of 
rotor-bladings I6j. Numerical solutions to subsonic flow by finite elements with 
self -adapting nodes given in Refs[l 6,1 9,22] show good agreement with experimental 
results. 

NOMENCLATURE: 



A, 

A, 

A2d 

A, 

Aa 

A t 

A s 

a 

C,W 

K,m 

M 


n 

P 



total area of boundary surfaces, 
inlet & outlet surfaces (Fig.l). 

periodic boundary surfaces (Fig.l): A 2u U A- d , A 2tJ - A^UA^ * 

A 2d -Ai d UAj d . 

free trailing vortex sheets, 
all solid boundary walls: A 3 = AciUA b - A3UA3*. 
hub- & casing annular walls: Aa^AiUAa . 
blade surfaces: A b -^A b UA b # . 

components of elementary area dA in r-, ¥>-, z- directions respectively, 

d A =d A V. 

shock surfaces, 
sound speed. 

absolute, relative flow velocity respectively. 

specific heat ratio C p /C v and mH/C-1) 1 
relative Mach number, 
outward normal unit vector. 

dimensionless pressure 
mass flux pA. 
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r 


r,<p,z and 

cylindnc coordinates fixed on rotor and the corresponding unit vectors. 
3ft, S rothalpy & entropy, 

s boundary position vector & its variation. 

u u-u/xr. 

dv elementary volume d v- rd<pdrdz. 

4> potential function 

U) angular speed of the rotor. 

A C ,A>A U demensionless forms of C,W,U respectively. 

P dimensionless density. 




Subscripts: 


annular and blade walls respectively. 
d,u downstream and upstream respectively, 

m circumferentially averaged value. 

n,T normal & tangential components. 

pr prescribed. 

radial, azimuthal and axial components respectively, 
p pressure blade surface, 

s suction blade surface. 

parameters on (A*, u , A* d ) and (A” u , A£ d ) respectively. 

- ,+ parameters just before & behind the interface or shock respectively. 

S uperscripts: 

o restricted variation 110]. 

*,** known and unknown portions of the boundary respectively. 

1. INTRODUCTION 


Nowadays the design of advanced turbomachinery would be impossible without 
using advanced aerodynamic theory and thereupeon based computational methods. 
During the last decade much progress has been made in this field and a detailed 
state-of-the-art review is given in Rer.[5], which reveals that with few exceptions, 
^•6* R^fs 11,2,11-17], most of work done to date, however, are concerned with the 
direct (analysis) problem, quasi-3-D flow model and mainly finite difference, finite 
volume and streamline curvature methods. Owing to the lack of exact (classical) 
variational principles (VPs) for rotor-flow finite element methods (FEM) used so 
far are exclusively based either on Galcrkin approach or on approximate VPs for 
the linearized problem. It is the great progress and the widespread and fruitful 
applications of the FEM in solid mechanics that motivated the present author in 
the mid 1970’s to start a systematic search for VPs in fluid mechanics in general 
I17,1SJ and in 3-D turbomachine flow theory in particular [1-4,8,11,12] with special 
emphasis on inverse and hybrid problems in order to provide both a new rigorous, 
sound theoretical foundation for FEM in computational aerodynamics of turbo- 
machinery and a number of novel rational versatile ways for new blade design or 
old blade modification. The hybrid problem is the one which, being a unification as 
well as generalization of the traditional direct and inverse problems, is capable of 
combining the merits of the two, while eliminating their shortcomings 11,2,11]. As a 
result, a lot of VP families have been established first for the direct problem 18] 
and the hybrid problem [11-13,16] of quasi-3-D cascade flow. In extending them to 
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fuly 3-D flow the major difficulties encountered were how to capture all possible 
unknown flow discontinuities, such as shock waves, free trailing vortex sheets and 
the unknown portions of the blade- and/or annular walls in inverse and hybrid 
problems. Subsequently, a series of VPs in terms of potential or stream functions 
for the direct problem of fully 3-D transonic potential and rotational flows with 
shocks in rotors have been developed in Refs [3,4 j, and furthermore, a unified 
variational theory of various hybrid problems for fully 3-D incompressible rotor 
flow has been presented iri Ref.ll] and extended to compressible flow in Ref .1.2], 
thereby a very powerful mathematical apparatus "the functional variation with 
variable domain" being used to full advantage for handling abovementioned flow 
diKcontinui ties. Successful numerical validations of such a theory have been 
carried out in Rcfs[l 6,1 9,22] by using a new finite element with self-adapting nodes. 

In the present paper, based on Refs [3,4], the unified variable-domain 
variational theory of hybrid problems for rotor-flow of Refs 11,2,16] is extended to 
fully 3-D transonic flow with shock waves in rotors of axial, radial- and mixed- 
flow f types. 

2 . BASIC : AERO T HERMODYNA Ml C_ EQU AT10NS 

Consider the fully 3-D subsonic and transonic potential, steady relative flow of 
an inviscid fluid past a rotating blading with constant angular speed u) (Fig.l). 

For such potential flows the nondimensional governing aerodynamic equations 
have the following form 13,4,7]: 

Continuity equation: 

v 4pA) J ( l ) 

r t dr dip dz > 
lrrotationality of the absolute flow: 

lip M (2) 

First law of thermodynamics: 

?+±iA 2 -Ai)-l (3) 

P 2m 

Homentropic equation: 


Eliminating p in Eq.(3) via Eq.(4) yields; 




Using Eqs. (2) & (3’), a full potential equation can be obtained from Eq.(l) 
[2,3,4,71: 


2M*M,s a # 'e'4, , 1+(M V +M„)' 3* 

- ■ — -=— ■ -2M r Mi- + = 

r 3z3v> 3r3z r 3r 


where 


2m 3r r 3 <p 


-AuYH^f -Ali 
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M ix > / ,■ l-'ni 

r — W 7 i 3 — / P y 

dr 

My,= W«> /a=( ^ ~ — A,.) / p' , 
r dip 

M.=W 7 /a-t— /P r ' r \ 

cz 

Mw=cjr/a- An ip' 

In the present paper, just as in Refs 11,2], the following three types of 3-D 
hybrid problems (H A XH A ), (H b \H A ) &. (H C XH A ) are studied in detail, where, as 
defined in Refs 11,2], e.g. (H c ;<H A ) denotes such a hybrid problem in which a hybrid 
problem of type C (H c ) is posed on the blade surface, while a hybrid problem of 
type A (H A ) is posed on the annular walls (Fig. 2a). In other words, the first symbol 

characterizes the problem type on the blade surface, while the second symbol 

that on the annular wall. As for the hybrid problems 1I A , 11^, , they arc defined 

in Table I for 2-D cascades 111,12,16], for the annular walls they are defined 
similarly 11,2]. 


Table I. Problem Classification (for the Blade Surface) 


i Types 

Given conditions | 

Geometric 

Aerodynamic 1 

H A 

Hb 

Part of airfoil form 

Pressure distribution along remaining part 

Airfoil thickness 
distribution 

i 

All cascade geometry 

Blade loading distribution (p P p s ) 

i 

" l i 

D ! 

Pressure distribution along suction surface 
Velocity difference distribution (yl s A P ) 
none 

i i 

i 

none 

L I 

Pressure distribution along airfoil contour j 


Of course, the abovementioned three 3-D hybrid-problem types arc only some 
typical ones taken as examples for consideration herein. Generally speaking, 3-D 
hybrid problems encompass a much wider variaty of types. They provides the 
designer with a scries oT novel design tools, which enable him to choose the most 
suited problem types or their combinations for meeting various practical design 
conditions at hand (e.g. aerodynamics, cooling, strength-vibrational and 
technological requiements etc). As pointed out in Ref .11], the three hybrid problems 
studied herein per se embrace also very comprehensive special cases, which, to a 
large extent, are capable of fulfilling various practical requirements of blade design 
and which can be made even much broader by posing different problem types on 
different portions of the blade (or annular) wall (Fig. 2b). 

3. VP F AMILY FOR THE (H A XH A )-P ROBLEM 

In the development to follow, starting from the VPs for the direct problem 
given in Refs [3,4] and employing the functional variation with variable domain 
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r 


[1,4], we present a unified variational theory of hybrid problems for fully 3-1) 
subsonic and transonic flows with shocks. 

Proceeding similarly to Ref •[! j leads to the following results without going into 
details. 

V P — 1« 1 he solution to the (H A • H A )-problem of 3-D sub- and transonic rotor-flow 

makes the functional J, stationary: tfj, -0, where <P y A*\ A /d , A. should be varied 
independent! y. 

I.+L+L aa , 


(6) 


where 


K, '‘llf pdv ‘I|j ! 

ivj I VI 

| |(qJp,^dA- II— 


1 — T" K|" )"+( rlf ) ?+ (|“ “ ~T“ 

2m dr r 3 <p 3^ r ■ 


3z 


3 <f' 


d v 


L- I j(qr i )pr^dA 

,A P 


A<P u )(pA n yd A 


L AA - 


| | (^W-dA+j | (?) 

.:*) .a**) 


) pr s dA . 


L 


With all unknown boundaries or interface A**> and treated by the method 
of functionnal variations with variable domain [1,4] the following set of 
stationarity conditions for J ; can be derived from 6J y 0: 

Euler’s eq.: Eq.(5) 

Natural boundary conditions (B.CJ: 
on A j : p A u- (q r*)pr , 

on A 2u : (pAnYHpAnY 9 , <P” - <P' A<P U , 

leading to the circumferential periodicity of all flow parameters, 
on A^ d ; {pAnY ipAn)^ ^, p’-p” , 

They arc just the interface conditions on the free trailing vortex 
sheets. 

on A s : Using <P_^<P+ as an essential (enforced) interface condition, we have 

O0/3r)_ ^-(3#/3 t) + , that is, the tangential velocity components at 
both sides of the shock are equal: 

Ut)-HAt)+ (7a) 

So we obtain the following natural interface conditions: 

(7b) 

(p//C+p/t£)_-(p//C+p.d£) + (7c) 

In addition, from Eq.(3) we can write 

(7d) 

Obviously, Eq-s (7a)-(7d) are none other than the well-known 
Rankine-Hugoniot shock relations [3,4]. 

on A^: pA n 0 . 

on AP: pA u = 0, and pm-~-(pm) pr . 

on A l*: pAn 0, and p^p pr . 

Thus, it has been shown that from this VP I actually the full potential 


i 


j 
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! equation (5) together with almost all boundary conditions for the 3-D (H A • H A )- 
problem can be derived naturally, and all unknown surface (e.g. shocks, free 
trailing vortex sheets and unknown walls) can be determined usin^ e.g. PKM. 

Applying a constraint-renu vmg transformation |17j, the above VP ! can be 
extended to the following generalized VP(GVP). 

JP fhc solution to the above 3-D (M A ■ H A )-problem makes the following 
functional J f/ stationary: with independent variations of <£>, A , p, /?, A - *, A 

and A^ d . 

T r (4>,A ,p,p,AV,A Sl A ?d M ;/ +l +k AA , (8) 


where 


p A . <P \( ' +A u +2A^A ¥ );+^P[ 1 ln(p : p K )}~ mpjdv 

In a way similar to the above one it can be shown that from dJ r/ 0 the 
following set of natural conditions results: 

Eulars equations: Eqs (l)-(4). 

Natural B.C.: All the same as those of VP L 

Via a constraint-recovering transformation [j'G, from GVP II a family of 
subgeneralized VPs can be derived, one of which is the foregoing VP I. 

4. VP I AMILY BOR/P HE (H b ;<H a )-PPOBL1{M 

in this case the B.C. on the annular walls remain the same as those 
(Ha v H a j-problem, while the B.C\ on the blade surface become: 

0) Blade thickness distribution given by 

‘Ar gf(r,z); 

(ii) blade-loading distribution given by 

Pr - V , x 
— £ — ^Kr(r,z) 

where g ¥ (r, z) and g f >(r,z) are prescribed functions. 

Proceeding in just the same way as in the foregoing section, we can establish 
the following VP family for the (H b XH A )-problem, which differs from that for the 
(h A XH A )-problem only in that the boundary integral term L AA now should be 
replaced by the following L BA : 


of the 

(9 A) 
(9B> 



[ BA__ 


| | ^ dA+| Jg, iTdA, 

(Aa*) (A t* } s 


( 10 ) 


imposing Eq.(9A) as an essential B.C.. In Eq.(10) the symbol (A to L stands for the 
suction blade surface. In this way we obtain the following VP family. 

VP 111: The solution to the (H B XH A )-problem of 3-D sub- and transonic rotor-flows 
makes the functions stationary: tf J m =0, thereby AT, A s , A h and should 
be varied independently. 


J r^r(^,Ar,A^, A,s,A2d)^E+L+L L 


( 11 ) 


It is easy to verify that from the same Euler’s equation and natural 

j B.C. set as those of VP I can be derived with the only exception that the natural 
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r B.C. on the blade surface A has now become Eq.(9B). 

GVP IV : The solution to t lie above 3-D (H b \H A )-probIem makes the following 

functional J ,, stationary: with independent variations of <P, A, p. p , A'*, A h , 


A- and A 


} n \<P,A f p f /> f A:* f A lf ,As-,A^)“l // +L+l. 


( 12 ) 


i 


From 6J lv 0 the same Euler's equations and natural B.C. set as those in GVP 11 
follow except only that the natural B.C. on the blade surface has become F.q. (9B). 

scy P ^Family: Bv means of a constraint-recovering transformation (171 a family of 
SGVPs can be derived from the GVP IV, including also the VP 111. 


5. VP FAMILY FOR THE (H c aH a )-PROBLEM 

Mere the B.C. on the annular walls still remain unchanged as before, while the 
B.C. on the blade surface arc now given as follows: 

(l) blade thickness distribution given by Eq.(9A): 

(li) pressure distribution along the suction blade surface: 

Ps- [p--(r,.:)] p , U3) 

To establish the VP family for the (H c X Ii A )-problen», wc proceed similarly as 
above. 11 turns out that this VP family differs from that for the (H A >.H A )-problcm 
only in that the boundary integral term l. AA should be replaced by the following 

L LA ; 


I. A - (^Vs-dA+f f(— s-dA, (14) 

lO 1A t>>s 

while the Eq.(9A) should be treated as an enforced B.C.. In Eq.(14) the superscript 
’o’ denotes that the 'restricted variation 1 * * * * 6 110] should be taken. Thus, we have: 

VP V : l or the (H c XH A )-prob!em 6*J^-0 with independent variations of 4 > , A*\ A D , A 
and A ?x1 holds, and 

J l ,(^,Ai*,A (E> ,A s ,A 2d )^I / +L+L CA f (15) 

VP VI: For the (H c XH A )-problem with independent variations of 4>, A , p, p, 

AJ\ A fc , As and A 2d holds, and 

,p,p,Ai , ,A b ,As,A 2 d) sas l// + I>L t A , (16) 

SGVP Pamilv : Applying the constraint-recovering transformation [17], we can 

derive a SGVP family for the (H c XH A )-problem from GVP VI, including also the VP 
V. 

It can be shown similarly as in previous sections that the Euler's equations 
and the naturall B.C. sets of the VP V, GVP VI and its derived SGVP family are 
the same as those of the VP I, GVP II and its derived SGVP family respectively, 
except that the natural B.C. on the blade surface now has become Eq.(13). 

6. SOME GENE RAL REMARKS 

1) It is easy to see that the traditional direct problem (DXD) [3,41 and inverse 


L 
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I” problem (1 >CI) are simply two special cases oT the (H A \H A )-problem, corresponding 
to A V 0 and A' -0 respectively. Accordingly, by setting A"-0 all VPs developed 
herein reduce to those presented previously in Kefs 13,4 J. 

2) IF, alternatively, the B.C. on the upstream periodic boundary A Zu (namely 
<P" <P'+-±<Pu) IS imposed as essential B.C'., the boundary integral terms on A 2u 
involved in L of all VPs should be dropped accordingly. 

3) An alternatve approach to handling free trailing vortex sheets is also 

possible by taking formally no variation of A 2d , though A 2d is unknown, but the 
interface conditions on A 2d (namely An^A'n^O, p^-p") are enforced as essential ones 
[ 21 ]. 

4) As stressed in Ref .11], sufficient attention should be paid to a rational 
choice of the position-variation 6s of the unknown boundaries A*/ and A^ d for 

facilitating the practical computation of (6s* dA) & 6J, (i I- VI). Some 

recommendations on this point are available in Ref .1 1 ] and, of course, also valid for 
the present case. 

For better shock-capturing a special finite element with self-adaptive build-in 
discontinuities is very promising and is now being under development. 

The numerical solutions to the problems I(H h +D)XD] & 1(H C +D)XD] obtained 
for subsonic flow by finite elements in Refs.[19,22] show good agreement with 
experimental results. 

8. CONCLUSIONS 

The unified theory of 3-1) hybrid problems of Refs [1,2] has been extended to 
transonic flow with shocks. This theory is pritnarly aimed at providing, firstly, a 
new rigorous theoretical basis of blade design for use in FEM and other direct 
variational methods (e.g. Ritz s method, Kantorovich's method) and, secondly, a wide 
variety of new rational versatile ways for new blade design and old blade 
modification. It also constitutes an important ingredient of the optimal design 
theory of 3-D rotor-bladings [6]. Based on the VPs for the direct problem of 3-D 
rotational flow [3,4], the present theory can be extended also to 3-D rotational 
flow. This will be presented in a companion paper. 
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